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Abstract 

We report a study of the modes DK^ and D*K^ where D^*^ decays to CP eigen- 

statcs. The data sample used contains 275 xlO^ BB events at the T(4S') resonance cohected by 
the Behe detector at the KEKB energy-asymmetric e+e^ coUider. The CP asymmetries obtained 
for DcpK are: Ai = 0.06 ± 0.14(stat) ± 0.05(sys), A2 = -0.12 ± 0.14(stat) ± 0.05(sys) and for 
D^^pK: At = -0.20 ± 0.22(stat) ± 0.04(sys), A^ = 0.13 ± 0.30(stat) ± 0.08(sys). 

PACS numbers: 14.40.Nci, 13.25.Hw, ll.30.Er, 12.15.Hh 
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Measurements of the decay rates of 5^ D^*'^K^ provide a theoretically clean method 
for extracting the Unitarity Triangle angle 03, an angle in the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix 0. Since both a and a can decay into the same CP 
eigenstate {Dcp, or Di for a CP-even state and D2 for a CP-odd state), the b ^ c and h ^ u 
processes shown in Fig. ^interfere in the DcpK^ decay channel. This interference 

may lead to direct CP violation. To measure D meson decays to CP eigenstates a large 
number of B meson decays is required since the branching fractions to these modes are of 
order 1%. To extract 03 using the GLW method the following observables sensitive to 
CP violation must be measured: the asymmetries 

^ B{B- ^ Di,2K-) - B{B+ ^ Di,2K+) 

" B{B- ^ Di,2K-) + B{B+ D,,2K+) ^ ^ 



2r sin S' sin i 



^3 



1 + + 2r cos 5' cos 03 
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and the double ratios 
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1^1,2 = -^po- = 1 + r + 2r cos 5' cos 03 (3) 

^' = \ ^ ^1 (A\ 
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The ratios P^^ ^ ^^d are defined as 



oDi,. ^ B{B- ^ Di,2K-) + g(P+ ^ 1^1,2^+) 
i3(P- ^ P'i,2vr-) + i3(P+ ^ P'i,27r+) 
^0 _ B{B- D'^K-) + B{B+ ^ 

i3(P- ^ + i3(P+ ^ D^Ti+) 

where r = |v4(P^ D^K^)/A{B^ D^K^)\ is the ratio of the magnitudes of the two 
tree diagrams shown in Fig. 5 is their strong-phase difference. The ratio r is given by the 
product of CKM factors and a color suppression factor, that characterizes the magnitude of 
CP asymmetry. The asymmetries and double ratios can be calculated for D* in a similar 
manner (notation A\2 and 7^12)- Here we have assumed that mixing and CP violation in 
the neutral D meson system can be neglected. 




u ^ u u ^ u 



FIG. 1: Feynman diagrams for B D^K and B — > D^K . 
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Previously, Belle j3| and BaBar |^ reported the observation of the decays B~ DiK~ 
and B~ D2K~ . BaBar also reported the observation of the decay B~ D^K^ . This 
paper reports more precise measurements of the DcpK~ channels, superseding our 

previous result, and a study of the decays B^ — > DIK~ and B~ — > D^K' with a data 
sample corresponding to 275 x 10^ BB pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [6]. 

mesons are reconstructed in the Cabibbo-favored modes (-D/) XII'- K^n^, CP-even 
modes (Di): K+R-, n-n+ and CP-odd modes (Da): fsT^TrO, K°u, K^g^). 

Neutral pions are reconstructed from pairs of photons, selected in the invariant mass 
range 118 MeV/c^ < M(77) < 150 MeV/c^, corresponding to a ±2.5(T window, where a is 
the vr" mass resolution. Each photon is required to have energy greater than 30 MeV in the 
laboratory frame, and also in this frame the pion candidate's momentum must exceed 0.5 
(0.1) GeV/c for the Kg7^^ i^s^) mode. The tt" candidates are kinematically constrained to 
the nominal mass. 

Each charged track not coming from a Kg candidate is required to be consistent with 
coming from the interaction point (IP). For each charged track, information from the ACC, 
TOF and CDC is used to identify the particle (PID) as either a pion or kaon via the K/ir 
likelihood ratio P{K/tt) = Lk/{Lk + -^^tt), where and are kaon and pion likelihoods. 
With the exception of the prompt kaon from the B meson decay ("fast track"), all kaon 
candidates must satisfy the PID requirement, P{K/7c) > 0.3. This requirement selects kaons 
with an efficiency of 92% and a pion misidentification rate of 18%. 

The candidates are formed from two oppositely charged pions with an invariant 
mass required to be within 8.5 MeV/c^ of the nominal mass (~ 3a). The meson is 
reconstructed from two oppositely charged kaons in a mass window \M{K~^ K~) —m^\ < 10 
MeV/c^. uj mesons are reconstructed from tt+tt^tt^ combinations in the mass window 0.757 
GeV/c^ < M(7r+7r"7r°) < 0.82 GeV/c^; a loose requirement of P{K/tt) < 0.9 is applied to 
the charged pions in the u candidate. 

For the D candidates, a 3a mass requirement is applied, where a is the D mass resolution 
which ranges from 5 to 12 MeV. D* candidates are reconstructed in the Dn^ decay channel 
depending on the decay channel. The mass difference between D* and D candidates is re- 
quired to be within 2.8 MeV/c^ (~ 3a) of the nominal value B meson candidates are 
formed by combining the D^*^ candidates with one charged track (denoted /i^). The signal 
is identified by two kinematic variables: the beam-constrained mass Mbc and the energy dif- 
ference AE calculated in the T(4S') center of mass (CM) frame, Mbc = \/-^beam ~ Ipd 
and AE = Ed + Eh ~ -Ebcam, where -Ebeam is the beam energy, pd and Eo are the momen- 
tum and energy of the candidate and ph and E^ are the momentum and energy of the 
K~ /tt~ candidate assuming the pion mass. With this definition, B~ — »• D^n~ events peak 
at AE = 0, while B~ D^K~ events peak around AE = —49 MeV. Signal candidates 
are selected with Mbc > 5.2 GeV/c^ and \ AE\ < 0.2 GeV. The experimental resolution for 
Mbc is approximately 3 MeV, dominated by the beam energy spread. The AE resolution is 
typically 10 MeV for all-charged-particle final states [Di modes). For final states with pho- 
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tons or neutral pions, the AE" resolution becomes broader and somewhat skewed to negative 
values. 

Event topology is used to distinguish BB events from continuum background. At the 
T(45'), the two B mesons are produced nearly at rest so these events tend to be spherical, 
whereas continuum events have a two-jet topology. We construct a Fisher discriminant 
of modified Fox- Wolfram moments called the Super- Fox- Wolfram (SFW) where the 
Fisher coefficients are optimized by maximizing the separation between BB events and 
continuum events. The angle in the CM frame between the B flight direction and the 
beam axis, cos 9b, is also used. These two independent variables {SFW and cos 9b) are 
combined to form a likelihood ratio: 71 = Lsig/(Lsig + Lcont), where Lgig and Lcont are 
defined as the product of SFW and cos 9b likelihood. The TZ requirement is optimized 
for each submode of DK and D*K using Ns/y/Ns + Nb, where Ns{Nb) is the expected 
number of signal (background) events in the signal region (the coefficients used for the Fisher 
discriminant are common to all the sub- modes). The expected number of simal events is 
obtained assuming the branching ratio given in the Review of Particle Physics [8] . Since the 
continuum background is negligible for the Kgcp mode, we do not apply an TZ requirement. 

For events with more than one candidate (l%-2% for all modes, except for KgU, ~10%), 
a single candidate is selected on the basis of a determined from the difference between 
the measured and nominal values of masses {D, D*, Kg, u, masses) and then the highest 
TZ value. 

Signal yields are obtained from fitting the AE distributions for 5.27 GeV/c^ < Mbc < 
5.29 GeV/c^. The PID for the fast vr or i^' is used to distinguish between D^*^tc and D^*^K 
modes (with a requirement P{K/tx) > 0.8 for D'^*^K, which selects kaons with an efficiency 
of 80% and a pion misidentification rate of 7%, and the remainder as D^*^7i). Signal peaks 
are fitted with double Gaussians. Shifts in the mean position and differences in resolution 
seen between Monte Carlo (MC) and data in Dn are used to correct the DK fits in data. The 
continuum background is modeled by a first order polynomial whose slope is obtained from 
the Mbc sideband (Mbc < 5.25 GeV/c^). Backgrounds from B meson decays are modeled 
by large MC samples using a smoothed histogram. When statistics are small, as is the case 
in D*K, shapes from D*7t are used directly. 

Backgrounds are studied using MC samples for known backgrounds and sidebands in 
data. A peaking background is found for B DK (where D tttt) coming from B Dn 
[D Kn), which is suppressed by making a 3a mass requirement on the Ktc invariant 
mass. 

For DK, in the K^K~ and tt+tt^ modes, clear peaks are seen in the D mass sideband, 
defined as 1.80 GeV/c^ < M{hh) < 1.83 GeV/c^ and 1.90 GeV/c^ < M{hh) < 1.93 GeV/c^, 
where h is a. charged kaon or pion. These peaks come from the B — > KKK and B — > Kirn 
modes, respectively. The yields obtained (63.5±7.5 events) are in agreement with the results 
of a dedicated study of these channels and allow an estimate of the peaking backgrounds 
for these modes. The sidebands are scaled (factor 0.5) and subtracted from the yields of Di, 
and hence are taken into account in the asymmetries and double ratios defined in Eqns. ^El 
Note that such effects are not seen in i? — D*K since the D* provides an extra constraint 
to reduce these backgrounds. 

Backgrounds in the D2 modes, K^uj and K^cj), need careful consideration because they 
can be modes of non-CP or with opposite CP (opposite asymmetry) to the mode considered. 
Possible backgrounds to K^uj include non-CP modes K^ttutt^ and K* p, and backgrounds to 
K^s'P include non-CP modes a^(980)/s:^, K^gKK, and opposite CP modes is:ga[](980) and 
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Kgfo{980). To determine contributions from these backgrounds, the data AE distributions 
for Dn modes are fitted in bins of or hehcity angle. The hehcity angle ^hei for (u) is 
defined as the angle between one of the kaons from (the normal to the uj decay plane) and 
D momentum in the (u) rest frame. The yields as a function of helicity angle are then 
fitted allowing for possible contributions from signal and either K^vrvrvr^ and K^KK. The 
fraction of signal is estimated to be 88.8 ± 8.4% for K'^u and 84.0 ± 12.5% for K^(f). When 
a helicity requirement (| cos6'hoi| > 0.4) is imposed to further reduce the backgrounds, these 
fractions become 92.4 ± 9.8% for K^u and 88.6 ± 11.1% for i^^0. 

The fitted AE distributions for positively and negatively charged B meson candidates 
are shown in Fig. |21 Table H] gives the corresponding yields and asymmetries with their 
statistical uncertainties. The asymmetries in the control samples (-D/) are consistent with 
zero, as expected. The modes of interest are DiK and D2K: the and B~ yields are used 
to calculate asymmetries after peaking background subtraction for Di. For D2 modes, the 
asymmetry is estimated mode by mode and the dilution factor due to the CP content of 
the background is taken into account, assuming no CP for K^uj background and opposite 
CP for Kg(f) background. 



TABLE I: Yields and asymmetries obtained for Dh and D*h modes. For D2 modes, the asymmetry 
is estimated mode by mode taking into account the CP content of the background. 





EB 


B+ 


B- 


A 


DfTT 


19266±150 


9677±103 


9521±102 


-0.008±0.008 


DiTT 


2163±56 


1049±38 


1124±37 


0.035±0.024 


D2Tr 


2168±61 






0.017±0.026 


DfK 


1131±41 


528±28 


603±29 


0.066±0.036 


DiK 


143.3±21.9 


70.2±14.7 


79.2±15.7 


0.060±0.144±0.046 


D2K 


149.5±19.0 






-0.117±0.141±0.049 


D*7r 


5434±101 


2756±59 


2678±59 


-0.014±0.015 


DItt 


662±37 


338±21 


322±21 


-0.021±0.045 


D*27, 


604±38 






-0.090±0.051 


D*K 


256±22 


140±16 


117±15 


-0.089±0.086 


DIK 


43.9±10.2 


27.3±7.4 


18.2±6.9 


-0.200±0.224±0.035 


D*2K 


32.7±10.0 






0.131±0.300±0.076 



The yields obtained for DIK and D^K (Fig.E)) are 43.9 ±10.2 and 32. 7± 10.0 respectively, 
which correspond to significances of 5.2a and 3.3cr {Kgir^ 2.9cr, KgU 0.9a, Kg(j) 1.4a) where 

the significance is calculated as y^— 2 ln(£o/>Cmax), where £max and £0 denote the maximum 
likelihood with the nominal signal yield and with signal yield fixed to 0, respectively. 

The sources of systematic errors for the double ratios come from the uncertainty in yield 
extraction, uncertainty in signal fractions for K^u and Kg(j) (1%) and the uncertainty in 
the contributions of peaking background from D sideband data. The uncertainty in yield 
extraction is estimated by varying the fitting parameters, such as the slope used for contin- 
uum or widths and means for Gaussians used for signals by ilex (6%-8%). The uncertainty 
due to peaking D sidebands is taken from the error on the estimated contribution: 6% for 
TZi. These errors are added in quadrature for d[*^ and D2*^ . 
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Systematic errors for A are from intrinsic detector charge asymmetry, measured from the 
control samples B — >■ -Djvr, (0.02), uncertainty in signal fraction for KgU and K^(f) and on 
the CP content assumption of the peaking background (0.01), yield extraction (0.02-0.04) 
and PID (0.01). 

The asymmetries for D12K, Ai and A2, are found to be: 

Ai = 0.06 ± 0.14(stat) ± 0.05(sys) 
A2 = -0.12 ± 0.14(stat) ± 0.05(sys). 

The double ratios are: 

7^l = 1.13 ± 0.16(stat) ± 0.08(sys) 
7^2 = 1.17 ± 0.14(stat) ± 0.14(sys). 



The asymmetries for Dl 2K are found to be: 

AI = -0.20 ± 0.22(stat) ± 0.04(sys) 
A*2 = 0.13 ± 0.30(stat) ± 0.08(sys), 

where the systematic errors are calculated in a similar way to the Dh case. The double 
ratios found are: 



7^* = 1.41 ± 0.25(stat) ± 0.06(sys) 
n*2 = 1.15 ± 0.31(stat) ± 0.12(sys). 



For each state such as B ^ DcpK four observables are measured (^1, A2, TZi, 7^2)- 
Since these are related by AiTZi = —A2'R-2 there are then only three independent observ- 
ables and there are also three physics quantities that should be extracted, 03, r^jx and Sbk- 
The measured asymmetries are consistent within errors with zero and also with the Stan- 
dard Model expectation. These measurements, while not sufficiently accurate to provide a 
measurement of 03, can be used to constrain 03 through a global fit jl2|, and the addition 
of modes such as D^pK can add further constraints 13l |. 



In summary, using 275 xlO^ BB events this paper reports results from the decays B^ — »• 
DcpK"^ and 5^ D^pK"^ where the CP eigenstates are those of the D meson. The 
partial rate asymmetries Ai^ are measured in both cases and are consistent with zero. The 
study oi B ^ DcpK is made with three times the statistics of previous studies. B D^K 
and B ■—>■ D2K are observed and their asymmetries and double ratios are also measured. 
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FIG. 2: AE distributions for (top left) 5+ DiK+, (top right) DiR-, (bottom left) 

D2K^, (bottom right) B~ D2K~ . Points with error bars are the data and the solid lines 
show the fit results. The components of the fit are the background from B meson decays (dotted 
line), the continuum background (dashed), the signal DK (left) and Dn (right). 
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